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Abstract 

The existence of traveling and standing waves is investigated for 
chains of coupled pendula with periodic boundary conditions. The 
results are proven by applying topological methods to subspaces of 
symmetric solutions. The main advantage of this approach comes 
from the fact that only properties of the linearized forces are required. 

This allows to cover a wide range of models such as Newton’s cradle, 
the Fermi-Pasta-Ulam lattice and the Toda lattice. 

Keywords: Newton’s cradle, coupled pendula, periodic waves, global 
bifurcation. MSC 34C25, 37G40, 47H11 

A chain of coupled pendula is a device of hanging limbs from an elastic 
rod. When this device is set in motion, each limb behaves like a pendu¬ 
lum that interacts with its neighbors by torsion forces. When the forces 
are approximated by Hooke’s law, the equations describing the pendula are 
equivalent to the discretized Sine-Gordon equations. 

Another device of interest is the Newton’s Cradle. This system consists 
of beads suspended by inelastic strings. In the absence of contact between 
beads, they have pendular motion; when beads collide with their neighbors, 
they repel each other with a Hertz’s type force. 

The paper aim is to study a model that includes the mentioned pendula 
among other situations of interest, such as the Klein-Gordon, Fermi-Pasta- 
Ulam and Toda lattices. 
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The movement of n coupled oscillators, qj(t) G M. for j = 1 with 

periodic boundary conditions, qj = qj +n , is described by equations 

- <?, = U\ qj ) + W\qj - qj-i) - W\q j+1 - Qj ), (1) 

where potentials U and W represent the dynamic and nonlinear interaction 
of the oscillators, respectively. 

The dynamic of a pendulum is governed by 

U(x) = 4(1 — cosx), uj 2 = cg/l ) (2) 

where g is the acceleration due to gravity and l is the length of the pendulum. 
The constant c represents the coupling strength after a rescaling, i.e. the 
normalized Hooke’s Law is given by 

W(x) = ix 2 , 

and the Hertz’s contact force by 

W(x) — - |*| 5/2 if x < 0, W(x) = 0 if x > 0. (3) 

5 

We assume that equations (JT]) have homogenous equilibria of the form 
qjit) = a for j — 1,.., n. 

These properties hold true in the coupled pendula and the Newton’s cradle 
when a = 0 and a = n. Under these considerations, the nonlinear equations 
(HD have periodic orbits arising from the homogenous equilibria. 

The study of periodic orbits in Hamiltonian systems goes back to Poincare. 
Lyapunov proves the nonlinear continuation of periodic orbits from normal 
modes of elliptic equilibria under non-resonant conditions. Later on, the 
Lyapunov center theorem was extended to consider the multiplicity of peri¬ 
odic orbits and the global properties of the families. Regarding multiplicity, 
the Weinstein-Moser theorem (1973) proves multiple periodic solutions of 
Hamiltonian systems with fixed energy, while the Fadell-Rabinowitz theorem 
(1978) considers fixed period. In [1J, Alexande-Yorke (1978) prove the global 
property of Lyapunov families. The proofs of these theorems make use of 
topological invariants that consider the 5'’-symmetry induced by time trans¬ 
lations. These theorems have been generalized to consider spatio-temporal 
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symmetries, for instance, see mm for fixed energy, mm for fixed period 
and mmm for global bifurcation. Other equivariant approaches can be 
found in j5[ |T2l 23] and references therein. 

In the context of equation (jT|) , in [IT] , local bifurcation of periodic solu¬ 
tions is proven in the case U — 0 and W" { 0) 7 ^ 0, and in pT3] , a Weinstein- 
Moser theorem is proven in the case U" = 1 and W"{ 0) 7 ^ 0. The present 
paper proves the bifurcation of periodic solutions using Brouwer degree in 
spaces of spatio-temporal symmetric functions, including a ^-symmetry in¬ 
duced by the reflection in time. In contrast with the results in m and [13] . 
this procedure has the advantage that allows to prove the global property. 

Theorem 0. Assume W" ( 0) 7 ^ 0, which is the case in the coupled pendula. 
For each k G [1, n/2) D N such that 

u k = \/U"(a) + ( 2 sin/c 7 r/n) 2 hF"( 0 ) > 0 (4) 

is non-resonant (Definition [7]), the homogenous equilibrium has three global 
branches of 27r/i'-periodic solutions; for k G {n/2, n}, only one branch exists. 
The frequencies v along the branches converge to v k as the solutions approach 
the equilibrium and the bifurcating branch is a continuum that either goes 
to infinity in Sobolev norm or period, or ends at other bifurcation point. 

In the coupled pendula, at the equilibrium a = 0, the non-resonant condi¬ 
tion of u k holds true except for a finite of parameters cu (Section 4.1). For the 
resonant parameters the theorem proves only the existence of the branches 
with the higher frequency Uj = lu k . 

When W"{ 0) = 0 all frequencies are resonant, 

v k \= \JU"[a ) for k G {1,..., n}. 

These resonances make impossible to obtain multiple periodic solutions by 
means of topological degree; instead, the existence of multiple standing waves 
is proven with the Fadcll-Rabinowitz theorem for odd potentials given in [ 6 ]. 

Theorem [Til Assume fF"(0) = 0 and U"(a ) > 0, which is the case in 
Newton’s cradle. The homogenous equilibrium has at least n/2 — 1 bifurca¬ 
tions of 27r/^-periodic solutions with symmetries (E]) and ([7]). The frequency 
v of the periodic solutions is arbitrarily close to n 0 , but the bifurcation does 
not necessarily form a local continuum. 

The periodic solutions have the symmetries of traveling and standing 
waves. This fact is proven exploiting the equivariance of equations ([T]) under 
the action of the group 

D n x 0(2), 
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where D n is composed by permutations of the oscillators and 0(2) shifts and 
reflects time; see Definition CD 

The symmetries presented in Section 3 are valid along the global branches, 
while the estimates shown in Section 2.1 are valid locally. We reproduce 
here the simple case n odd and k — 1. The symmetries for k — 1 have 
been analyzed previously in m and [T2]. The present paper completes the 
classification of the symmetries for all k’s. 

Symmetries and estimates. Let 

( = 27T /n, qj(t) = a + Xj(ut), 

where Xj(t ) is 27r-periodic and v is the frequency. The branch of traveling 
waves has symmetries and local estimates, 

Xj{t) = Xn-j(-t) = x j+1 (t - C), (5) 

Xj(t ) = r cos (t + j() + 0(r 2 ), 

where r is a parameterization of the branch and 0(r 2 ) is a 27r-periodic func¬ 
tion. The branch of standing waves has symmetries and local estimates, 

Xj{t)=x n - j {t)=x j {-t), ( 6 ) 

Xj(t) = rcos(jC) cost + 0(r 2 ), 

and the other, 


Xj (t + 7 t ) = X n -j{t) = Xj(-t), (7) 

Xj(t) = rs'm(j() sin t + 0(r 2 ). 

Traveling waves for Newton’s cradle have been estimated asymptotically 
in [TE] . and for beads in na and [22]. I 11 these papers, traveling waves 
are constructed by means of a reduction to a single equation with delay. 
This procedure is commonly used in many problems; see |T9] and the ref¬ 
erences therein. However, the reduction to one equation cannot be used to 
prove existence of standing waves, and then, one of the achievements of the 
present paper is the construction of them for the Newton’s cradle. Neither 
the Weinstein-Moser theorem in [[13] is applicable to the Newtons’s cradle be¬ 
cause it assumes non-resonant conditions over hh"(0). Through Proposition 
E] a Weinstein-Moser theorem can be proven for the Newtons’s cradle. 

In the case of beads, W"(0) = 0 and U = 0, the topological approach 
used to establish the existence of standing waves cannot be used due to the 
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fact that Vfr — 0 for k — 1, n. In H3 is shown that standing waves exist 
even for homogenous potential W with = 0. Therefore, further work is 
necessary to investigate the existence of standing waves in the beads problem. 

In Section 1, we set the bifurcation problem and make a global reduction 
to a finite number of Fourier components. In Section 2.1, we prove the 
global bifurcation in the case W"{Q) ^ 0. In Section 2.2, we prove existence 
of standing waves in the case W"{ 0) = 0 and U"(a) > 0. In Section 3, 
we describe the symmetries. In Section 4, we apply the theorems to the 
coupled pendula and the FPU and Toda lattices. In Section 4.3, we present 
a comment about the existence of standing waves in the case of a homogenous 
potential W with = 0. 

1 Setting up the problem 

Let q = (qi, ... ,q n ) and 

n 

r(?) = £[t/te) + w%,(8) 

3 = 1 

Equations (P) can be expressed in vectorial form as 

-? = w(?). 

Hereafter, we assume that the potential V is twice differentiable and 

a = (a,...., a) 

is an equilibrium, VV(a) = 0. This is equivalent to assume that U\a) = 0 
and W'{ 0) = 0. 

Using the change of variables q(t) = a + x(ut), the system of equations 
become 

— u 2 x = VU (a + x). 

Let H? h (W : ) be the Sobolev space of 27r-periodic functions. We define the 
operator / from if| jr (M n ) in Ll n (M. n ) as 

/(i;y) = -A-W( a + x). (9) 

Since a is an equilibrium, then /(0; u) — 0 for all v. Therefore, the branches 
of 27r/z/-periodic solutions emanating from the equilibrium a correspond to 
zeros of f(x] v) bifurcating from (0, u Q ). 
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Definition 1 Let D n be the subgroup of permutations generated by 

((j) = j + 1 , n{j) = n — j modulus n. 

Let p : D n x 0(2) —> GL(L\ n ) be the homomorphism generated by 

p(l)(xi,:.,x n ) = (x 7 (i),...,x 7 ( n )), ( 10 ) 

for 7 G D n and for ip, R G 0(2), 

p(ip)x(t) = x(t + ip), p(R)x(t) = x(—t). (11) 

Then p defines a D n x 0(2)-representation o/L 27 r (M n ) and induces the left 
action in L^fMT) given by p( 7 , x) = p(^)x . 

Since V(x) is invariant by the action of D n , then VV(x) is D n -equivariant. 
Thus, the operator f(x) is D n -equivariant. Given that the equations are 
autonomous and reversible in time, then f(x) is D n x O(2)-equivariant. 

1.1 Lyapunov-Schmidt reduction 

In the case W"( 0) 7 ^ 0, we will prove existence of periodic solutions using 
a global Lyapunov-Schmidt reduction and Brouwer degree. Although an 
application of Leray-Schauder degree can provide similar results without 
reductions, we prefer this approach because the reduction is required in the 
Newton’s cradle anyway. The idea of the global Lyapunov-Schmidt reduction 
is taken from [15]. 

The Fourier expansion of x G L\ % and the projection P are defined as 
x(t) = ^r^xie llT and Px = xie M . 

IGl'Zj |Z|<Z 0 

Let Xx and x 2 be the components of x, given by 

X! = Px, x 2 — (I — P)x. 

The components of / are given by f± = Pf and / 2 = (I — P)f. 

We realize the global Lyapunov-Schmidt reduction in the set Q p x A e , 
where 

Ll p = {x G : \\x \\ H 2 < p}, A £ = {v > e}. 

27T 
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Proposition 2 There is a Iq such that x 2 (xx,zz) is the only solution of 
/ 2 (xx + x 2 ,z/) = 0 in Tip x A e . Thus /(xx + x 2 , is) — 0 if and only if 
</>(xi, is) = 0, where 

= /i(xx + x 2 (xx,i/);i/) (12) 

= -is 2 d tt x i - PX7V (a + xx + x 2 (xx, z/)). 

Furthermore, the reduced map 0(xx,zz) is _D n x 0(2)-equivariant, where the 
action of ip,R E 0(2) in i/ze l-th Fourier component is given by 

p((p)xi = e dlfi xi, p{R)xi = xi. 

Proof. If we find a positive constant a such that 

||<9 X2 / 2 (xx + x 2 )y 2 || L 2 >a\\y 2 \\ H 2 

Z7T 2 7T 

for all (x, is) E Tl p x A e , the global implicit function theorem due to Hadamard 
(Theorem 5.1.5 in [3]) implies existence of a unique function x 2 (x 1; z/) such 
that / 2 (xi + x 2 (z/, xx), is) = 0. Using the uniqueness of x 2 (xx,z/), it can be 
proven that /i(xx + x 2 (xx, is), is) is D n x 0(2)-equivariant. 

Let y 2 G (/ — P)H'i n . Since V G C 2 (M n ) and ||x|| c o < cp for all x G Tl p , 
then ||0 2 U(a + x)|| c o < Cp for all x G Tl p . Therefore, 

\\(I-P)D 2 V(a + x)y 2 \\ L2 <Cp\\y 2 \\ L% < C(p/i 0 2 ) fe|U . 

11 11 ^2-n 2-7T 27r 

Using the previous estimate and \\is 2 d u y 2\\ L 2 > H 2 / 2 llixr 2 f° r u e we 

27r 27r 

have 

11 <9x 2 /2 2/2 11 L 2 > ||z/ 2 a tt 2 / 2 || 2 - ||(J-P)U> 2 U(a + x)|/ 2 || 2 >o ||?/ 2 ||„2 , 

where a = e 2 — C(p/ll). We conclude that a > 0 if Zo > Cy/p/e. Note that 
the number of Fourier components Iq goes to infinity as e —* 0 and p —* 00 . 


Using Taylor’s expansion of the map / 2 (xx + x 2 ), we obtain the estimate 



||x 2 (xx,z/)||< c||xx|| 2 

(13) 

for x close to 0. Thus, the linearization of the reduced map at (0, is) 

is 


</>'(0; z/)x 1 = -is 2 d tt x 1 - P> 2 U(a)xx = M(lis)x h 

(14) 


\i\<i 0 


where 

M(lis) = (lis) 2 I-D 2 V( a). 
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1.2 Irreducible representations 

In this section we identify the irreducible representations of 0(2) x D n . In 
the Z-th Fourier component, the action of 0(2) is 

p(p)xi = e hip xi, p(R)xi = xi. 

The Fourier components are subrepresentations of the group 0(2); then, we 
need to find the irreducible representations of xi G C n under the action of 
D n . 

Let ( = 2tt/ n and e k G C n be 

e fc = n~ 1/2 (e 1{ikC \ e 2{iK \ ..., e n(ifcC) ). (15) 

The vectors e k for k = 1 are orthonormal and their direct sum is the 
whole space C n . Therefore, we can expand xi = Yhk&z ri Xk F w h ere ^n = 
{1, n}, and 

x(t) = ^ x k je k e llt . 

(k,l)€ Z„xZ 

Proposition 3 For k = nj 2, n the action of the group D n x 0(2) in x k ,\ G C 
is given by 

p( (,¥>) x k,i = ±e l(p x ktl , p(K)x k> i = x ki i, p(R)x kt i = x kfl , (16) 

with negative sign fork = n/2. Fork G [l,n/2)flN, the action in (x kj i,x n _ ki i) G 
C 2 is given by 

p(c, <p){x k ,i, Xn-k,i) = e lip (e lkC x k ^ e~* K x n _ M ) 7 (17) 

P(^)(x k) li X n — k \) (Xn— k ,liX k \), 

p(n)(x k i 1 x n — k) \) (x n —k,iiiX k i). 

Proof. The actions of ( and k in e k are 

p(C)e k = n~ 1/2 (e 2{ik<: \ e 3{ik ° ,..., e n( * fcC) , e (ifc °) = e ik<: e k , 
p( K )e k = ( e 0-m^ e {n-2 )tfcc ?e 2*K e iifcC ? e ™*C) = e n _ fc . 

Moreover, the action of k G 0(2) is 

p(k)e k z = n- x ' 2 (e- iK z, e~ 2iki z, ..., e^-^z, e“ nifc ^) = e n _ fc z. 

The result follows. ■ 

By the previous proposition, the subspaces generated by e n , e n /2 and 
e k © e n - k for k G [l,n/2) fl N are irreducible representations. 


1.3 Isotropy groups 

In this section, we identify the maximal isotropy groups of the irreducible 
representations presented in Proposition [3j 

The action of D n x 0(2 ) in the subspaces generated by e n and e n / 2 are 
given by where the minus sign corresponds to the representation gen¬ 
erated by e n / 2 - The subspace x n ^ E K has isotropy group 

Tn = {( C,0),(«,0),(0,k)), (18) 

and the subspace x n /2,i E M, 

T n /2 = ((C,^), («,0), (0,«)). (19) 

Both isotropy groups T k for k = n, n/2 have fixed point spaces of dimension 
one. 

Definition 4 Let h be the maximum common divisor of k and n, 

1=^ a ndn = ^. 
h h 

The cases k E [l,n/2) D N are analyzed in the appendix. We reproduce 
here the main results. 

Proposition 5 For k E [l,n/2) n N, the representation e k © e n - k has three 
maximal isotropy groups with fixed point spaces of dimension one. An isotropy 
group is 

T k = (((, ~k(), (k, k), (n(, 0)) . 

For h odd, the other isotropy groups are 

Sk = ((«, 0), (0, R), (■ n(, 0)), S k = ((k, t r), (0, n tk), «, 0)). 

The case h even is given in the appendix. 

These isotropy groups are relevant because the reduced map 

f H (x-v) : Fix(H) -E Fix(ff) 

is well defined and under non-resonant conditions, the linearization Df H ( 0, v) 
has a simple eigenvalue in the kernel. 
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1.4 Linearization 

The Hessian of V is 


D 2 V( a) = U"(a)I + W"{G)A, 

where A = (a t j)™ J=] is the matrix defined by = 2 if i — j, a UJ = —1 if 
\i — j\ — 1 modulus n, and a M - = 0 otherwise. 

Proposition 6 T/ie matrix M{y) = i^ 2 / — T> 2 P(a) is diagonal in the basis 

j 

x= 22 M(v)x = 22 X k (v)x k e k , (20) 

where the eigenvalues are 

X k (u) = u 2 - U"(a) - (2 sin k(,/2) 2 W"(0). (21) 

Proof. Since the j-th coordinate of e k is n _1 / 2 e®-^, then the j-th coordinate 
of Ae k is 

(2 - (e _ifcC + e ikC ))n- 1/2 e iK = 4sm 2 (k(/2)n~ 1/2 e iK . 

We conclude that Ae k = 4:sm 2 (k(/2)e k and 

M(y)e k = [u 2 I - U"{a)]e k - W"(0)Ae k = X k (u)e k . 


The fact that X n - k = X k for k e [l,n/2) D N is consequence of Schur’s 
lemma. That is, since e k © e n - k is an irreducible representation, Schur’s 
lemma implies that M(v) = XI in e k © e n - k . 

In the basis x k ,i, we have 

</>'(0; ^)xi = 22 ^k{lv)x K ie k e llt . (22) 

\l\<l 0 k^^n 

2 Main results: Bifurcation theorems 

We proceed to prove two bifurcation theorems. 
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2.1 Coupled pendula 

If W"{f)) ^ 0, we avoid resonant frequencies by assuming that D 2 V(a) is 
invertible, which holds if U"(a) and W"{ff) are positive. 

Definition 7 We say that the frequency 

v k = y/U"(a) + (2 sin kn/n) 2 W"(0) (23) 

is non-resonant if lv k 7 ^ Vj for j G (k,n/ 2] D N and l > 2. 

Theorem 8 Assume W"{fS) 7 ^ 0 and D 2 V( a) is invertible. For each k G 
[ 0 ,n/ 2 ] fl N such that u k > 0 is non-resonant, the equilibrium a has three 
global bifurcations of 2tt/ v-periodic solutions emanating from v = v k with 
isotropy groups T k , S k , and S k . 

Proof. The linear map 0'(O; v) has Z-th Fourier block M(lv) with eigenvalues 

\j{lv) = {Ivf - v*, 

for j G [0, n/2] fl N. 

For l = 0, the matrix M(0) = —D 2 V( a) is invertible by hypothesis. For 
l > 2 , since W"{fS) ^ 0, frequencies Uj are increasing in j. Thus lu k > Uj for 
j G [0, k\ fl N and, by hypothesis, lv k ^ Uj for j G (fc, n/2] fl N. Therefore, 
matrices M(lv k ) are invertible for l > 2. Given that M(lv ) is continuous in 
v, the blocks M(lv ) are invertible for v close to v k . 

For 1 = 1, since u k is increasing in k G [0,n/2] flN, then v\ — u 2 ^ 0 
for j 7 ^ k. Therefore, the linearization <f>'( 0; v) is invertible for v close to v k , 
except for the block that corresponds to the representation (x kil ,x n _ k ^) G 
e k © e n — k , 

A k (u)I : C 2 —> C 2 . (24) 

Set H equal to T k , S k or S k . Since the group H has fixed point spaces 
of real dimension equal to one in e k © e n - k , the restriction D(p H ( 0; v k ) has a 
simple eigenvalue crossing zero in (l24lh Using Brouwer degree as in Theorem 
14 in [9] or [15], we conclude existence of a local bifurcation in the fixed point 
space of H. 

The global property follows from assuming that the branch is contained 
in the set Qjf x A e for some e and p, unless it is an unbounded continuum 
set with period or Sobolev norm going to infinite. Applying Brouwer degree 
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to the reduction in x A e , we conclude that the sum of the local degrees 
at the bifurcation points is zero, as in Theorem 5.2 in |15j or Theorem 15 in 

i- ■ 


Remark 9 In the case U — 0, the Hessian D 2 V( a) has a zero-eigenvalue 
corresponding to the conserved quantity Y^j=i Qj = 0- If this is the only zero- 
eigenvalue, we can extended the previous theorem using the restriction of f 
to the subspace 

n 

x = {xe l? 27T : = 0 }. 

3 = 1 

Also, analogous theorems can be proven in the case of D n -equivariant long 
range interactions, W(qi, ...,q n ). 

Under non-resonant considerations, the local reduction can be realized on 
the 1-th Fourier component 


E ilt 
xie . 

I*l=i 


Let r be a parameterization of the amplitude of the local branch, where 
the frequency v is a functions of r. The branch with isotropy group T k has 
eigenvalues corresponding to the coordinates (x k ,i,x n _ ki i) = (r, 0) (see the 
appendix). We conclude that the projected component x x can be estimated 

by _ 

xi (t) = r(e lt e k + e it e k ) + 0(r 2 ). 

Moreover, since e lt e k + e lt e k = 2 cos (t + jkCf) and x 2 (x x ; v) = 0(r 2 ), then 

Xj(t) = 2 r cos (t + jk() + 0(r 2 ), (25) 

where 0[r 2 ) is a 27r-periodic function of order r 2 . 

Similarly, the coordinates for the isotropy group S k are {x kt i,x n _ kt i) = 
(r, r), then 


x x (t) = r(e lt e k + e u e k + e lt e n _ k + e it e n _ fc ) + 0{r 2 ). 

We conclude 

Xj(t) = 4r cos jk( cos t + 0(r 2 ). (26) 
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For the group Sk and n odd, the coordinates are (xk,i,x n -k,i) — (n ~ r ), 
then 

x j(t) = — 4r smjkfsmt + 0(r 2 ). (27) 

For n even, the coordinates are (xfc,i, Xn-fc, 1 ) = (r, re*^), then 

x,-(t) = —4r sin(jfcC — C/2) sin(t + C/2) + 0(r 2 ). (28) 

Note that for n odd, since sin(27rj/n) = 0 for j e {n, n/2} D N and 
cos(27 ij/n) 7 ^ 0 , standing waves Sk and Sk have nodes of different orders, r 
and r 2 respectively. 

2.2 Newtons’s cradle 

For Newtons’s cradle, W is given by fl3J) and U(x) = o; 2 (l — cosx). Since 
Vk = oj for k = 1 ,..,n, the previous theorem cannot provide the existence 
of many periodic solutions. In this case, the application of Weinstein-Moser 
theorem guaranties the existence of at least n periodic solutions in each 
constant energy surface, 

H(q,p) = ^ ||p || 2 + V(g). (29) 

However, these solutions may agree with the n traveling waves found in jT 6 j , 
and the symmetries have to be considered in order to obtain new solutions. 

Definition 10 Let S = S 1 and S be the group generated by (k, tt) (0, ttR) if 
n is odd, and 0 ) and (O,C«0 if n is even. 

We will prove existence of standing waves using the fixed point spaces of 
the groups S and S, which are the isotropy groups Si and Si without the 
generator 

Since Uk = cu for all k, the 1-Fourier component is non-resonant with 
other Fourier components. Thus the local reduction of Section 1.1 can be 
realized on the 1-th Fourier component, 

xi = xie lt + xie~ lt , 

and the local bifurcation map is defined by (^(xp^) : C n —> C n for (x,is) 
close to ( 0 , oj). 
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Since the operator f{x) is the gradient of 

F(x) = J \d t x\ 2 -V(x)^j dx, 

we can conclude (see Section 1.9 in [15]) that 0(xi) is the gradient of the 
reduced potential 

$(xi) = F( Xl + x 2 (xi)). 

That is, $'(x!)/i = (</>(x 1 ),h ) i 2 . 

Since Lyapunov-Schmidt reductions preserve equivariant properties, the 
potential $(xi,v) is D n x 0(2)-invariant. Setting H equal to S or S, this 
implies that the restriction Q H : Fix(H) — y R is well dehned and the gradient 

is 

V<f> H = (f) H : Fix(TT) Fix(if). 

Furthermore, the element n G 0(2) is contained in the Weyl group of H 
and acts multiplying by —1 the fixed point space of H (see appendix), then 
<h H (xi) = <h^(— Xi). 

Using a = 0, U(0) = 0, O 2 U(0) = lu 2 I and x 2 (xx;z/) = 0(|xx| 2 ), we 
estimate 

<E>(xi, v) = 27 t(v 2 - u 2 ) |xi|“ + o(|aq| 2 ). (30) 

Therefore, the Fadell-Rabinowitz theorem in [6] implies that the odd po¬ 
tential $ H (xi,i/) has at least dim Fix(H) branches of critical points. Since 
dimFix(if) > n/ 2 — 1 for H equal to S and S (see the appendix), the fol¬ 
lowing theorem holds. 

Theorem 11 Set H equal to S or S. If (xi, v) = (0,o;) is an isolated point 
of the potential u), then f(x, u) has at least d_ zeros for u < c 0 and 

d + for u > co in the fixed point space of H. The zeros converge to (0,cu) as 
v —y co and 

gL + d + > n /2 — 1. 


Therefore, there are at least n /2 — 1 bifurcations of periodic solutions 
with symmetries S and S emanating from the homogeneous equilibria. The 
information about the symmetries of these nj 2 — 1 bifurcations can be im¬ 
proved if one considers similar procedures in the fixed point spaces of Sk and 
S k . 
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Another proof of the previous theorem can be given with Z 2 -equivariant 
Conley index. The idea in [3j consists on using the invariant property of 
the Conley index and the change in the unstable set of the gradient flow 
generated by , which has dimension dimFix(if) for v < oj and 0 for 
v > OJ. 

Using P rop osi t ion ITT a ml Theorem 9.9 in [3], we can prove that the energy 
surface H~ 1 (e) for small £ contains at least nj 2 — 1 solutions with isotropy 
groups S and S (Weinstein-Moser theorem). These kinds of solutions are 
called brake orbits in [3]. A Weinstein-Moser theorem for coupled pendula 
was proven in [13] using results of [18] - 

3 Description of symmetries 

The isotropy groups T k , S k and S k have the generator (n£,0). Functions 
fixed by (n(, 0 ) satisfy 

Xj(t) = X j+n (t), 

where h is the maximum common divisor of k and n, k = k/h and n = n/h. 

Therefore, the oscillators in these solutions form a wave of length n that is 
repeated h times along the pendula. In the following discussion, we describe 
only the wave of length n, 

Oi, ...,x n ). 


3.1 Traveling waves 

The group T k has generators and solutions with isotropy 

group Tfc have symmetries 

Xj(t) = Xn-j(-t) = x j+ i (t - k(27i/n)), (31) 

For k = 1, two consecutive oscillators have a phase shift of 2n/n, 

Xj(t) = Xj + i(t — 2tt /n). 

For k 7 ^ 1, the solutions are just permutations of the case k — 1. 
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3.2 Standing waves of the first kind 

Standing waves have different behavior depending on the parity of n and 
nj 2. We present three cases for each kind of standing waves. Given that 
cases k 1 are permutations of k = 1, we present only the case k — 1. 
Hereafter, we denote 

X*(t) = Xj(t + 7r). 


n odd 


Since the isotropy group Sk has generators (k, 0) and (0 ,k), functions fixed 
by Sk have symmetries 

Xj(t) = Xn-j(t ) = Xj(-t). 

Setting n = 2m + 1, the wave of length n is 


(xQj X \, ..., X m , X m , ..., Xi), 


where Xj are even functions for j = 0, ...,m. 


(32) 


n/2 odd 

The isotropy group Sk has generators (k, 0), ,tt) and (0, R), solutions in 
the fixed point space of Sk satisfy 

Xj{t) = X n -j{t) = x* j+n/2 (t ) = Xj(-t). 

Setting n = 4m + 2, the wave of length n is 


('Y */-. O'* ^ rf rp ^ rp ^ rp * rp ^ rp ^ /y* O'* . \ 

\ U/ 0') •••? •••? *^0’ ^1’ **•’ ^mi j 


(33) 


where x,- are even functions for j = 0, m. 


n/2 even 

Setting n = 4m + 4, we conclude that the wave of length n is 


l 'Em. 


* * * 
rp rp rp 
l *^1 1 -hfll 'T 1 






(34) 


where Xj are even functions for j = 0, ...,m and x m is 7r-periodic. 

Observe that, although solutions (|33li and (f34li have the same isotropy 
group, the solutions are qualitatively different at the oscillator x m . 
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3.3 Standing waves of the second kind 

n odd 

The isotropy group Sk has generators and (0,7 tk). Functions fixed by 

S k have symmetries 


Xj(t) = Xfi—j (t + 7r) = 

Setting fit = 2m + 1, the wave of length n is 

(xo, X\, ..., x m , x m ,. 

where Xj{- + vr/2) are even functions for j 
periodic. 


Xj{TT-t). 

(35) 

= 1 and x 0 (t) is even n- 


n/2 odd 

Since the isotropy group S k has generators (tv£, 0 ), (|£, 7 t) and ( 0 , ^/c), solu¬ 
tions in the fixed point space of S k satisfy 

Xj(t) = Xn-u+1 )(t) = x*n/ 2 +j(t) = x j {2n/n - t). 

Setting n = 4m + 2, the wave of length n is 

(x 0 , x u x m = x* m ,x*, x* Q , x* Q , x* m = x m , xi, x 0 ), (36) 

where Xj(7i/n + •) are even for j = 0, ...,m and x m is 7 r-periodic. 

n/2 even 

For n = 4m + 4, we conclude that solutions with isotropy group S k satisfy 

(xo, ■ ■ • 5 x m , x m iX q, Xq,.. x m , x m ,...j Xq), (37) 

where Xj(n/n + •) are even for j = 0,m. 

Solutions (Effij) and (IH7D are different close to the oscillator x m . 


4 Applications 

In most applications the potential W is convex. If U is also convex, the 
frequencies v k are always positive; if U is concave, the frequencies u k are 
positive for k G [k 0 ,n/2] fl N. 
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4.1 Coupled pendula 

The coupled pendula via torsion springs are governed by the dynamics of 
U(x) =uj 2 (1 — cosx) andVh(x) = x 2 /2. Since U'(0) = oj 2 sin x and VF"(0) = 
1, the homogenous equilibria are 0 = (0,0) and n = ( 7 r,7r). 

Given that U"( 0) = c o 2 and U"(tt) = —co 2 , the condition that gives the 
existence of bifurcation is 


v k = \J ±o; 2 + (2 sin kir/n ) 2 > 0, 


where the minus sign correspond to n. For the 0 equilibrium, the frequencies 
v k are always positive. For the 7r equilibrium, despite the fact that the un¬ 
coupled system is unstable, the coupled system has 27r/z/fc-periodic solutions 
near the equilibrium for k G [ko,n/2] D N, where 


ko ~ (u/tt) arcsin(o;/2). 


In Theorem [ 8 j the non-resonant condition of v k is equivalent to oj 2 7 ^ 
±Ui(j) for integers / > 1 and j > k, where 


ui(j) ■ = 


(2 sin kir/n) 2 — (2 sin jn/n) 2 /! 2 
1-1 /l 2 


(38) 


For the 0 equilibrium, the non-resonant condition c o 2 7 ^ uji (j ) holds true 
except for a finite number of resonant parameters c oi(j), and for n, except for 
countable number of parameters —c oi(j) — > 2 sin knfn as l —> 00. 

The same statements are true for the discrete Klein-Gordon equation 
with potentials U = u> 2 x 2 + x 3 and W{x) = x 2 /2 in [20], and for the bistable 
potential U{x) = u> 2 ( 1 — x 2 ) 2 /4 for equilibria with a = 0 and a = ±1. 


4.2 FPU and Toda latices 

In this section, we consider oscillators with U(x) = 0 and nonlinear interac¬ 
tions 

OO 

W (x) = x 2 /2 + PkX k - 

k =3 

The FPU lattice corresponds to W{x) = x 2 /2 + (3x 3 /3 and Toda lattice to 
W(x) = e~ x + x — 1. 
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Since IU'(0) = 1, then 0 is an homogeneous equilibrium. Therefore, for 
each k 6 [1, n/2] D N such that 

Uk = 2 sin kir/n 

is non-resonant, the equilibrium has three global bifurcating branches of 
periodic solutions. The non-resonant condition Uj ^ is equivalent to 
ui(j) 0. 

Actually, in [2Tj is proven that resonances (uj = lu k ) and higher order 
resonances do exist. In these cases, Theorem [S] only proves the existence of 
the bifurcations with the higher frequency Vj — lv k - hi [IT] a two dimensional 
family of superposed traveling waves is constructed for resonances {uj = Ivk) 
but in the context of infinite FPU lattices. 

4.3 Homogeneous lattices 

In this section we discuss the existence of standing waves when u k = 0 for all 
fc’s. We reproduce the case U = 0 and W(x) — | \x | 5 ^ 2 presented in [IT] , 
Equations ([Tj) have standing waves of the form q 3 (t) = a,jq(t ) when 

- q = W'(q ) (39) 

and 

— a,j = W'(dj + 1 — a,j) — W'{dj — Oj_i). (40) 

Since W (x) is convex, equation (139|) has only periodic solutions. 

Let bj = W'(a,j — a 3 -\ ) be the momentum of a 3 . Equation (140]) is equiv¬ 
alent to 

—cij = bj .|_i — bj, cij — cij—i = bj \bj\ / . 

Let 0 : M 2 M 2 be 

4>(a, b) = (a + b \b\~ 1/3 ,b — (a + b |6| _1/3 ) j . (41) 

The orbits (a 3 ) 3e z of 0 are solutions of (14H . 

In [TS] is shown that the map 0 is conservative and has many invari¬ 
ant orbits near (0,0). These invariant orbits correspond to periodic and 
quasiperiodic (in space) standing waves. These arguments exhibit existence 
of standing waves even in the case that topological methods cannot be ap¬ 
plied. Therefore, further investigations are required to establish existence of 
standing waves for the beads problem studied in [IT)] and [22] . 
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Appendix 

The action of D n x 0(2 ) in the irreducible representation (z \, z 2 ) G ei© e n _i 
is given by 


p( C, ¥>)(zi, 22 ) = e vp (e' < z 1 , e~ zC z 2 ), (42) 

/»(«)(*!, * 2 ) = (%,2l), 

P(k)(2i,2 2 ) = 

Proposition 12 T/ie representation has three maximal isotropy groups 
with fixed point spaces of real dimension equal to one. The isotropy groups 
are: 


Parity 

Orbit points 

Generators 

Isotropy group 

n = any 

(r,0) 

(C,-0>(«>«) 

T\ 

n = odd 

(r,r) 

(r, -r) 

(«,0), (0,k) 

(K,7r), (0,7m) 

Si 

Si 

n = even 

(r,r) 

(r, re *^) 

(«,0), (0,«), (7T,7r) 

(«c,0), (0,C«), (7r,7r) 

Si 

Si 


Proof. We need to analyze isotropy groups of orbit types. By applying k, we 
can assume that z\ 0, unless both coordinates are zero. Using the action 
of S 1 , we can assume that (zi,z 2 ) = (r, pe 1 ^) with r > 0. Moreover, by the 
action 

p(K-K)(r,pe i *) = (r,pe<*-^), 

we have that if 6 (—(, £], and by the action of (k, r) that 6 [0, £]. There¬ 
fore, we need to analyze only isotropy groups of 

(zi,z 2 ) = (r,pe i for if e [0, C] ■ 

Actually, for n odd, the element 2( generate the group Z n and we can take 

[0,C/2]. 

The elements (l(, (p), (nl(, gf) G D n x S 1 act as 

p(K,p,)(r,pe*) = (re'^\pe < »-‘^), 
p(nl(,ip)(r,pe‘*) = (pe i{ *~ li+r \ re‘ K+v) ). 
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The first coordinate r is fixed by (ZC, ip) if p = — ZC, and by (nl(, ip) if p = 
l( — 0 mod 2tt. These elements act as 

p(i C ,-K )(r, P e*) = (r,pe‘»- 2 '«), (43) 

p(K/C, - ip)(r,pe'*) = (p, re'i 2K ~' t ' 1 ). (44) 

Thus, we need to find out when these elements fix the second coordinate. 

For p = 0, elements (ZC, —ZC) always fix (r, 0). For p 0 0, elements 
(ZC, —IQ fix ( r 5 pe l Q if -0 — 2/C = 0 or /C = vr (mod 27t). That is, the element 
(7T, 7T) is in the isotropy group of (r, pe l Q when n is even. 

For p 0 r, elements (k/C, ZC — 0) never fix these points. For p — r, 
elements — 0) fix (r, rewhen 2/C — 0 = 0 (mod 27 t) or 0 = ZC 

(mod 7r) . Then (k0, 0) is in the isotropy group depending on the parity of n. 
If n is odd, the element (k, 0) fixes 0 = 0, and (k, it) fixes the point 0 = 7r. 
If n is even, the element (k, 0) fixes 0 = 0, and (kQ 0) fixes the point 0 = C- 
We conclude that the orbit points in the table are fixed by the generators. 
Moreover, the fixed point spaces of the isotropy groups are the set of orbit 
points (in the table) for rgK. Therefore, in the irreducible representation, 
the fixed point spaces have real dimension equal to one. ■ 

The action of D n x 0(2) in the irreducible representation e*, © e n -k is 
given in Proposition |3l This action is similar to (1421) . except that C £ D n 
acts as 

d(C)0i,-2) = (e lfcC £i,e" lfcC z 2 ). 

Proposition 13 Fork G [l,n/2)DN, the representation ek®e n -k has three 
isotropy groups with fixed point spaces of real dimension equal to one. The 
isotropy groups are: 


Parity 

Orbit P 

Generators 

Isotropy G 

n = any 

(r,0) 

(C,-fcC),(K, «),(wC,o) 

T k 

n = odd 

(r, r) 

(r, -r) 

(«, 0), (0,k), (nC,0) 

(«, 7r), (0, ttR) , (nC,0) 

s k 

S k 

n = even 

(r,r) 

(r, re im Q 

(«,0), (0,k), (fmC,vr), (nQ 0) 
(kjtiC, 0), (0, m (k), (fmC, tt), (n(, 0) 

S k 

s k 


where m G N is the modular inverse of k (mk = 1 mod n). Furthermore, the 
Weyl groups of these isotropy groups contain the element tt G S 1 that acts 
multiplying by — 1 the fixed point spaces. 
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Proof. The action of £ G D n can be expressed as 

p(C)(zi,z 2 ) = (e iH2ir/n) z 1 ,e~ lH27T / n) z 2 ). 

The factor k acts as a permutation because k and h are coprime numbers. 
Since mk = 1 modulus h, then 

p(m C) = (e l ( 2ir / n \ e _ * (2,r/n) ), p = ~I, p(n C) = I. 

The results in the table can be obtained using the previous proposition with 
m( = ‘lix jn instead of (. Furthermore, the element n G 0(2) leaves invariant 
the fixed point spaces and acts as — 1. Then i r is contained in the Weyl 
groups. ■ 

The 1-th Fourier mode X\ G C n is the direct sum of the irreducible rep¬ 
resentations ei, e n /2 and e*, © e n -k for k G [l,n/2) D N. Let 

Fix(iL) = {x\ e C n : p(y)xi = 11,7 6 id}. 

Proposition 14 For n odd, 

dim K Fix(S') = n/2 + 1 / 2 , dim® Fix(S') = n/2 — 1 / 2 , (45) 

and for n even, 

dim R Fix(S') = n /2 + 1 , dim K Fix(S') = n/ 2 . (46) 

Proof. For n odd, the group S is generated by («, 0) and (0, R) and S by 
(k, 7 r) and (0, ttk). For k G [1, n/2)flN, the point (^ 1 , ^ 2 ) G efc®e n _fc is fixed by 
S if z\ = z 2 G M and by S if z\ = — z 2 G R. Then, the fixed point space of the 
groups S and S have dimension one in ek © e n -k ■ We conclude that the total 
dimension of the fixed point space for the representations k G [1, n/2) fl N is 
(n — l)/2. For k = n, the point x n< i G C is fixed by S if x n< i G R, and by S 
if x Ut i = 0 . 

For n even, the group S is generated by («</, 0) and (0, C^) • F° r k G 
[1, n/2) fl N, we have the actions 

p(0,C«)(^i,^ 2 ) = (e* c ^ 2 ,e* c zi), 
p( K C; 0)(zi, z 2 ) = {e~ lK z 2 ,e lK zi). 
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Thus, the point {z\, z 2 ) G e*, © e n _fc is fixed by S if z 1 = z 2 G M, and by 5 if 
Zi = e lt ^z 2 and Z\ = e^ k+1 ^z 2 . Both conditions define subspaces of dimension 
equal to one. Therefore, the total dimension of the fixed point space, for 
k G [l,n/2) fl N, is n/2 — 1. For k = n/2,n, the point Xk,i G C is fixed by S 
if Xk t i G M. Moreover, the representation k = n is fixed by S if x Ut i = re 
■ 
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